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ABSTRACT 

Single  tensile  overloads  were  applied  to  4340  steel  specimens  which 
were  heat  treated  to  give  120  and  220  ksi  yield  strength  levels.  The 
influence  of  yield  strength  level  on  the  number  of  ncnsteady  state  crack 
growth  cycles  subsequent  to  the  application  of  r.  ! O'J  percent  overload 
was  noted  to  be  substantial.  The  number  of  nonsteady  state  cycles  for 
the  120  ksi  yield  strength  steel  was  approximately  an  order  of  magnitude 
greater  than  that  of  the  higher  strength  steel.  A  retardation  model  was 
developed  using  a  residual  stress  intensity  factor  concept  similar  to 
that  prcocscd  by  Willenborg  et.  al.  The  model  was  found  to  predict  to 
within  10  percent  the  number  of  nonsteady  state  crack  growth  cycles 
required  to  move  a  crack  from  the  pre-overload  position  to  a  subsequent 
position,  one  overload  induced  plane  stress  plastic  zone  radius  ahead  of 
the  pre-overload  position.  The  model  indicates  that  the  reason  for 
substantial  increases  in  nonsteady  state  crack  growth  cycles  observed  for 
the  low  strength  steel  is  due  to  a  corresponding  increase  in  the  overload 
affected  zone  size. 
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SECTION  I 

INTRODUCTION 

Variable  amplitude  load  spectra  when  applies  to  cracked  structures 
result  in  observations  of  crack  movement  which  are  at  variance  with 
simple  no-load-interaction  crack  growl!'  prediction  techniques.  One 
recognized  facet  of  the  crack  growth  behavior  generated  under  spectrum 
loading  is  that  substantially  slower  crack  growth  rates  are  observed 
when  small  amplitude  load  cycles  are  preceded  by  a  load  level  with  a 
higher  maximum. 

Several  crack  growth  retardation  models  have  been  developed  to 
predict  this  high-low  load  interaction  effect.  Wheeler  (Reference  1) 
provided  a  model  which  identified  the  concept  of  a  load  interaction  zone. 
He  derived  this  model  assuming  that  the  ratio  of  interaction  zones  is 
directly  related  to  the  observed  decrease  in  crack  growth  rates. 

Willenborg  et  al  (Reference  2)  used  the  Wheeler  formulation  of  load 
interaction  zones  to  derive  a  model  which  was  based  on  the  difference 
between  load  interaction  zone  sizes.  Both  of  the  above  mentioned  models 
were  reformulated  by  Gallagher  (Reference  3)  to  help  identify  the  physical 
implications  as  well  as  to  provide  identifiable  modeling  forms  for  others 
who  are  directing  their  attention  to  a  residual  stress  intensity  factor 
concept. 

The  use  of  a  residual  stress  intensity  factor  concept  suggests  that 
the  local  or  effective  stress  intensity  factor  which  is  sensed  by  the 
propagating  crack  be  calculated  from 

Keff  -  k"  -  Kr  (1) 

where  K  is  the  stress  Intensity  factor  associated  with  remotely  applied 
loads.  The  residual  stress  intensity  factor  (K^)  would  in  general  be 
considered  a  function  of  many  interacting  variables,  e.g.: 
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where 


KD  =  f  (K0L  ,  K0L  ,  Nn.  ,  K™  ,  R°°,o  ,fT ,  ...) 
R  'max  min  OL  max  ys 


(2) 


=  maximum  K  generated  by  the  overload 
max  3  J 


=  minimum  K*  generated  by  the  overload 


number  of  overloads  applied 

mavimum  generated  by  the  loads  following  the  overload 


R  =  stress  ratio  of  remotely  applied  loads 

U.  GO 

N  =  number  of  cycles  of 
o  =  yield  strength  level 
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SECTION  II 

TEST  METHODS  AND  MATERIALS 

1.  MATERIALS  AND  SPECIMEN  GEOMETRY 

The  nominal  chemical  composition  of  the  one-half  inch  thick  4340 
steel  specimens  which  were  employed  in  this  investigation  is  given  in 
Table  I  and  the  heat  treatment  schedules  used  to  produce  the  220  ksi  and 
120  ksi  yield  strength  levels  are  given  in  Table  II. 

The  geometry  of  the  specimens  is  shown  in  Figure  1  and  the  corres¬ 
ponding  stress  intensity  factor  coefficient  is  given  in  Figure  2  from 
Reference  4. 

2.  TEST  TECHNIQUES 

A  constant  stress  intensity  factor  condition  was  held  throughout  the 

OL  / — ; — — 

tests  such  that  an  overload  stress  intensity  factor  1^  of  40  ksi  v  in 
was  followed  by  repeated  cyclic  load  applications  corresponding  to  a 
remotely  calculated  stress  intensity  factor  varying  between  2  and  20 
ksi  v'  in.  The  condition:  ^  *  40  ksi  /"7n,  =  20  ksi  /in  , 

K^in  *  2  ksi  /in' was  applied  repeatedly  to  specimens  heat  treated  to 
the  two  yield  strength  levels  while  observations  of  crack  length  movement 
were  made.  Prior  to  a  reapplication  of  an  overload,  it  was  determined 
that  the  cracking  rates  associated  with  the  K^ax  *  20  ksi  /TrT condition 
were  steady  stote  and  no  longer  appeared  to  be  influenced  by  the 
previously  applied  overload. 

A  schedule  of  loads  for  any  0.025  inch  increment  on  the  crack  path, 
as  required  by  Figure  2  to  achieve  the  desired  constant  stress  intensity 
conditions,  was  developed  prior  to  test  start  up  and  then  utilized  during 
the  test.  It  was  normally  issible  to  define  the  surface  crack  position 
to  within  a  O.OOOo  inch  tolerance  since  crack  length  measurements  were 
made  using  a  20X  Gaertner  microscope  mounted  on  a  1.0  inch  Aerotech 

translation  stage  which  was  remotely  controlled  to  move  in  0.0001  in:h 
increments. 
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TABLE  I 


NOMINAL  CHEMICAL  COMPOSITION  OF  4340  STEEL 


Si 

Ni 

Cr 

Mo 

Fe 

0.28 

TABLE  II 

HEAT  TREATMENT  SCHEDULE  FOR  4340  STEEL  SPECIMENS 
FOR  220  KSI  YIELD  STRENGTH 

1.  Austenitize  1600°F  for  1  hour  in  salt  bath 

2.  Oil  Quench 

3.  Temper  1  hour  at  600°F  in  air 

4.  Water  cool 
Repeat  steps  3  and  4 

6.  Surface  grind,  remove  0.020  in.  from  both  sides  of  the  specimen 

FOR  120  iGMTaDTrRENGTH-” 


1.  Repeat  steps  1  thru  6  (220  k si  yield  strength) 

2.  Anneal  at  1 550° F  for  1  hour  in  air 

3.  Slow  cool  specimen  in  sand  from  1550^  to  ambient  temperature 
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All  Dimensions  In  Inches 
Tolerance  1 0.005 


Figure  1.  Geometry  of  Single  Edge  Notch  (SEN)  Specimen 
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Loads  were  applied  to  the  specimen  with  a  mechanical  test  systerr 
(hydraulically  actuated,  servo-controlled)  which  had  a  load  cell  mounted 
in  series  with  the  specimen.  The  frequency  of  applied  loading  was  10  Hz 
for  the  constant  amplitude  (sinusoidal  wave)  K^ax  =  20  ksi  /Tn, 

K^in  =  2  ksi  /TrT  cyclic  condition  and  approximately  0,05  Hz  (manually 
applied)  for  the  overload  applications.  The  high  frequency  loading  and 
a  desiccated  air  environment  were  employed  to  limit  the  influence  of  any 
moisture  induced  environmental  action  (Reference  4). 


7 


AFFDL-TR-74-27 


SECTION  III 

RESULTS  AND  DISCUSSION 


1.  CRACK  GROWTH  BEHAVIOR 

Detailed  in  Figure  3  is  a  schematic  showing  the  parameters  that  have 
been  utilized  recently  (References  5,  6,  and  7)  to  characterize  the  over¬ 
load  affected  crack  propagation  behavior.  The  nonsteady  state  (overload 
arrowed)  region  is  defined  in  terms  of  cycles  (N*)  and  distance  (a*), 
termed  the  overload  affected  crack  length,  over  which  the  reducing  effect 
is  observed. 


The  observed  crack  growth  behavior  of  a  crack  which  has  been  subjected 
a  single  100  percent  overload  (k^L  =  40  ksi  /  in  ,  =  20  ksi  /"Tn) 

is  presented  in  Figures  4  and  5  for  the  120  and  220  ksi  yield  strength 
levels,  respectively.  The  data  are  presented  in  terms  of  an  incremental 
movement  of  the  crack  as  measured  from  its  position  immediately  prior 
to  the  overload  application.  The  data  shown  in  Figures  4  and  5  were 
utilized  to  find  N*  and  a*  for  each  overload  application.  Table  III 
contains  a  listing  of  the  observed  N*  and  a*  parameters  and  the  steady 
state  crack  growth  rates  measured  outside  the  delay  region. 

The  average  a*  and  N*  parameters  calculated  using  the  data  listed 
in  Table  III  directly  indicate  the  influence  of  yield  strength  on  over¬ 
load  delay  behavior.  For  the  120  ksi  yield  strength  steel,  the  average 
overload  affected  crack  length  is  16.6  mils  which  is  traversed  in 
approximately  33,000  cycles;  the  corresponding  values  for  the  220  kri 
yield  strength  steel  are  4.6  mils  and  4100  cycles.  Or  the  / s  of  this 
information,  one  would  have  to  rank  the  120  ksi  y it  ■  •t.i.j'h  steel 

superior  to  the  220  ksi  yield  strength  steel  iur  spectrum  load  usage  due 
to  the  significant  increase  in  the  number  of  reduced  growth  rate  cycles 
observed  in  this  steel.  Furthermore,  the  average  steady  state  crack 
growth  rates  for  the  120  ksi  yield  strength  steel  are  approximately  a 
factor  of  2.5  lower  than  those  associated  with  the  220  ksi  yield  strength 
level . 
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Figure  5 
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TABLE  III 

OVERLOAD  AFFECTED  ZONE  PARAMETERS 


jOl 

max  _  0 

CO 

K 

max 

120  KSI 

YIELD  STRENGTH 

220  KSI 

YIELD  STRENGTH 

0VFRL0AD 

a* 

mils 

N* 

Aa 

?  ss 

10  D  in/cycle 

mils 

N* 

Aa  1 

-6  ss 

10~°  in/cycle 

1 

13 

25,500 

1.50 

6 

3800 

4.71 

2 

16 

34,500 

2.06 

4 

5600 

4.69 

3 

19 

32,500 

1.22 

5 

4000 

4.14 

4 

15 

30,000 

0.87 

4 

4000 

3.75 

5 

17 

35,750 

1.44 

5 

4900 

3.70 

6 

(7 

33,500 

1.83 

4 

2600 

3.12 

7 

19 

37,000 

1.31 

4 

3700 

3.05 

3 

17 

36,000 

1.16 

5 

4400 

2.85 

9 

- 

- 

- 

4 

3900 

2.66 

AVERAGE 

16.6 

33,100 

1.43 

4.56 

_ m _ 

is  the  growth  rate  following  overload  affected  growth  region, 
ss 
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2.  OVERLOAD  DELAY  AFFECTED  CRACK  LENGTH 

As  discussed  by  Gallagher  (Reference  3),  both  Wheeler  (Reference  1) 
and  Willenborg  et  al  (Reference  2)  assumed  that  the  overload  affected  crack 
length  would  be  a  fraction  of  the  overload  established  load  interaction 
zone  (ZqL).  Wheeler  assumed  that  Zql  could  be  approximated  by  the  over¬ 
load  created  plane  strain  plastic  zone  radius  while  Willenborg  et  al 
utilized  the  overload  plane  stress  plastic  zone  radius  to  characterize 
Zq^.  The  corresponding  estimates  of  the  overload  affected  crack  length 
used  by  Wheeler  and  Willenborg  et  al  are 


(3) 


(4) 


respectively, 


Recently,  several  investigators  (References  5,  6,  and  7)  have 
employed  variations  of  the  following  relationship  (r  is  the  plastic 
zone  diameter  for  the  applied  overload)  to  show  good  correlation  with 
measured  overload  affected  crack  lengths 


*  r 


1 


o  '  2X+6CPXT 


(5) 


where  Prcbst  and  Hi  11  berry  (Reference  7)  assumed  K+  *  and  both 

Von  Euw  et  al  (Reference  5)  and  Trebules  e*  al  (Reference  6)  assumed 

K+  *  -  K01:  ,  and  X  represents  the  fraction  of  the  cross  section 

max  min 

exhibiting  shear  lips.  From  the  fracture  surfaces,  X  was  determined 
to  be  near  zero  for  both  yield  strength  levels  for  this  investigation: 
Equation  5  reduces  to  „  ,  /  „+  \  2 


% 


1 

y? 


for  the  case  when  X  is  zero 


© 


(6) 
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Table  IV  provides  a  direct  comparison  between  the  average  overload 
affected  crack  lengths  (a*)  determined  from  the  data  given  in  Table  III 
and  the  parameters  calculated  by  Equation  3,  4,  and  6  with  K+  =  K^. 

It  appears  on  the  basis  of  this  comparison  that  the  Willenborg  et  al 
plane  stress  assumption  provides  the  more  viable  approximation  to  the 
affected  crack  length. 


TABLE  IV 


COMPARISON  OF  MEASURED  AND  CALCULATED  AFFECTED  ZONES 


Jys 

(KSI) 

a* 

avg 

(in.) 

r 

0 

(in.) 

r* 

J'OL 

(in.) 

%L 

(in.) 

120 

0.0166 

0.0118 

.00625 

0.0177 

22G 

0.00456 

0.00350 

.00186 

0.00525 

Note  that  the  inverse  ratio  of  the  two  corresponding  yield  strength 
levels  squared  (220/120)  closely  predicts  the  ratio  of  the  average 
measured  affected  crack  lengths  (16.6/4.6).  It  is  expected  that  future 
observations  will  further  substantiate  this  influence  of  yield  strength 
(in  a  given  material)  on  the  size  of  the  affected  crack  length. 

3.  WILLENBORG  ET  AL  MODEL  PREDICTED  CYCLES 


One  test  for  any  crack  growth  predictive  model  is  its  ability  to 
predict  the  number  of  cycles  of  applied  load  to  achieve  a  given  crack 
length  change.  The  Willenborg  et  al  model  as  detailed  by  Gallagher 
(Reference  3)  takes  a  form  identical  to  that  given  by  Equation  1  where 
the  residual  stress  intensity  factor  KR  is  defined  by 


yPl 

max 


1/2 


max 
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where  Aa  is  the  crack  growth  into  the  overload  load  interaction  zone 
(Zq^).  Following  an  overload,  Aa  is  approximately  zero  and  the  reducing 
stress  intensity  factor  is  maximum.  For  the  test  conditions  considered 
in  this  investigation,  Equation  7  gives  an  initial  of  20  ksi  /in. 
Equation  1  therefore  predicts  that  the  maximum  effective  stress  intensity 
factor  is  zero  at  the  crack  tip. 

The  maximum  effective  stress  intensity  factor  must  be  greater  than 

the  threshold  (maximum)  stress  intensity  factor  (approximately  6  ksi  v^n 

for  steels  when  R  =  0)  in  order  that  fatigue  crack  growth  occurs.  The 

assumption  expressed  by  Equation  7  must  be  in  error  since  Table  III  shows 

that  a  finite  number  of  cycles  was  required  to  propagate  the  crack  through 

the  overload  affected  crack  length.  Others  (References  5  through  11) 

have  also  observed  that  cracks  grow  after  single  overload  applications 
0L  00 

where  K  /K  *  2.0,  the  condition  for  which  the  Willenborg  et  al 
max  max 

model  predicts  zero  crack  tip  (effective)  stress  intensity  factors. 

4.  ALTERNATE  CRACK  GROWTH  DELAY  MODEL 

A  more  general  approach  to  developing  the  residual  stress  intensity 
factor  would  suggest  that  KR  may  be  proportional  to  the  Willenborg  et  al 
residual  stress  intensity  factor  expressed  by  Equation  7. 

Kr  =  4  Kj  (8) 


Consider  the  appropriate  test  conditions  required  to  evaluate  the 
proportional i ty  factor  $.  One  boundary  condition  might  be  the  "shut-off" 
overload  to  maximum  load  ratio  that  produces  no  crack  growth. 

This  would  be  established  using  techniques  similar  to  those  employed  by 
Probst  and  Hi  11  berry  (Reference  7).  Assume  that  the  shut-off  overload 
level  develops  a  local  stress  intensity  factor  condition  such  that  no 
growth  is  induced.  Since  the  fatigue  threshold  stress  intensity  factor 
maximum  (K  )  is  approximately  constant  for  negative  stress  ratios 

"loXyu 

"  eff 

(Reference  12),  set  the  maximum  local  stress  Intensit"  factor  (K  ) 

*  max 
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equal  to  K  for  zero  tension  loading  (R  =  0).  Using  Equations  7  and 

mat  X  *r  U 


,eff 


TH 


8,  K  can  be  expressed  as 
max 


„eff 

"max 


,  I  yO L  n  Aa  J/2  i/00 

max  '  v  1  max  ^  Zq^'  max 


(9) 


Immediately  following  the  shut-off  overload  K^x  which  produces  no  growth, 

:.a  =  0  and  K  =  K  and  Equation  9  can  be  solved  for  o  subjected  to 
max  max-m 

these  assu  iptions: 


K  '  -  K 

max  max 


TH 


'ni  to 
Cx  -  Kmax 


(10) 


Probst  and  Hill  berry  (Reference  7)  presented  results  in  which 

increasing  higher  overloads  eventual ly produced  a  condition  of  no  crack 

growth  for  five  different  \C  levels  in  2024-T3  Aluminum.  They  found 

max 

that  the  shut-off  overload  ratio  was  given  by 


K°l 

max 


ax 


2.3 


OD 


for  each  k"  level  studied, 
max 
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5.  CRACK  GROWTH  PREDICTIONS  WITH  NEW  MODEL 


In  order  to  develop  c  estimate  of  the  number  nonsteady  state 
cycles  (N*)  associated  with  the  overload  affected  zone  the  crack  grov/th 
rates  that  result  from  each  application  of  loading  must  be  integrated 
over  the  overload  affected  zone.  A  Paris  power  law  equation  (Reference 
13) 


(12) 


was  chosen  as  an  adequate  representation  of  the  crack  growth  rate  data 
for  stress  ratios  less  than  0.1  (assuming  that  tension  only  stress 
intensity  factors  establish  .IK). 

The  steady  state  crack  growth  equations  used  for  predicting  N*  were 


1-376  x  lO'10  (13) 

for  the  120  ks 1  yield  strength  steel  and 

^  =  3.46  x  10'10  £K (14) 


for  the  220  ksi  yield  strength  steel.  The  exponent  3.2  was  determined 
from  Independent  tests  on  these  two  steels  by  others  (References  4  and 
14)  and  the  pre-exponential  constant  was  selected  so  that  Equations  13 
and  14  would  intersect  the  average  steady  state  crack  growth  rates  at 
K  a  18  ksi  /Tn  given  in  fable  III. 


The  effective  stress  intensity  ranges  used  in  Equations  13  and  14 

eff  eff  eff 

were  obtained  either  by  subtracting  K  ^  from  K^  when  K  ^  was  greater 


„eff 


eff 


than  zero  or  by  the  value  of  K^  when  K^  calculated  using  Equations 


1  and  9  was  less  than  zero. 
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The  parameter  ;  defined  by  Equation  10  was  estimated  on  the  basis  of 
two  assumptions:  (1)  the  shut-off  overload  ratio  for  the  steels  investi¬ 
gated  is  close  to  that  defined  by  Equation  11  and  (2)  K  for  steels 

maxTj_| 

is  approximately  6  ks i  /TrT  for  a  stress  ratio  of  zero.  Equation  1C 
subjected  to  these  two  assumptions  becomes 


1  -  6/K 

O 


max 


(lb) 


A  numerical  integration  approach  was  employed  to  estimate  the  nu-ber 

of  applied  cycles  required  to  propagate  the  crack  from  its  preoverload 

position  to  a  position  one  plane  stress  (overload)  yield  zone  radius 

(Equatior  4)  aheac  of  the  preoverload  position.  One-tenth  of  the  plane 

stress  yield  zone  radius  was  thought  to  be  an  adequate  increment  in 

the  integration  scheme.  Table  V  lists  the  calculations  corresponding 

to  the  development  of  the  predictions  based  on  the  new  model  defined  by 

Equations  1,  7,  8,  15,  and  13  or  14.  Table  V!  summarizes  the  results 

of  a  *,  N*.  r  ,  N**  (average  number  of  N  measured  to  achieve  r  1 
y0l  y0L 

and  the  corresponding  predictions  associated  with  both  the  Willenborg 

et  al  and  new  models.  As  can  be  discerned  from  Table  VI,  the  new  model 

predicts  (to  within  10  percent)  the  number  of  cycles  observed  to  achieve 

r^  .  The  new  model  is  also  noted  to  be  conservative  since  steady  state 

rates  are  predicted  in  a  shorter  number  of  cycles  than  was  observed. 


6.  PREDICTION  OF  CRACK  GROWTH  RATES  USING  NEW  MODEL 

T oe  crack  growth  data  presented  in  Figure  4  were  used  to  obtain  the 
cyclic  crack  growth  rates  for  the  120  ksi  yield  strength  4340  steel 
listed  in  Table  VII.  The  fatigue  crack  growth  rate  data  for  the  first 
five  singularly  applied  overloads  are  presented  graphically  in  Figure  6. 
While  the  magnitudes  of  the  crack  growth  rates  associated  with  each  curve 
plotted  are  noted  to  be  similar,  the  e  appears  a  tendency  for  the  shape 
of  the  trough  in  the  curve  to  become  flatter  and  longer  with  each  succeed¬ 
ing  overload.  This  observation  is  reinforced  when  Table  VII  is  used  to 
compare  the  growth  data  for  the  first  and  second  overload  to  that  tor 
the  seventh  and  eighth  overload.  The  reason  for  this  change  in  crack 
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TABLE  VII  (Continued) 


Based  on  Linear  Extrapolation 
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growth  behavior  with  each  successive  overload  is  due  to  the  thru-the- 
tnickness  crack  front  curvature  changes  which  will  be  discussed  below. 

Table  VIII  provides  a  summary  of  the  crack  growth  rate  behavior  for 
the  first  five  singularly  applied  overloads  as  a  function  of  the  crack 
movement  subsequent  to  the  applied  overload.  These  data  were  then  used 
to  obtain  the  average  growth  rates  as  a  function  of  position  for  the  first 
four  and  first  five  overloads.  Table  V  lists  the  crack  growth  rates 
computed  using  the  model  as  a  function  of  incremental  position  in  the 
overload  affected  zone.  Figure  7  presents  the  average  growth  rate  for 
the  first  four  and  first  five  overloads  (see  Table  VIII)  and  the  calcu¬ 
lated  crack  growth  rates  for  the  new  model. 

Figure  7  shows  that  the  observed  behavior  is  substantially  different 
from  the  model  prediction  in  that  the  observed  growth  rates  reduce  grad¬ 
ually  to  their  lowest  level  while  the  model  predicts  an  abrupt  change 
in  growth  behavior  followed  by  a  monotoniccl ly  increasing  growth  rate. 

It  can  be  noted  that  the  model's  monotonically  increasing  crack  growth 
rate  parallels  the  observed  average  growth  rate  behavior  as  it  again 
increases  to  its  steady  state  level.  If  the  steady-state  crack  growth 

region  at  the  end  of  the  r  (i.e.  4.7  mils)  is  shifted  to  a  position 

y0L 

prior  to  the  beginning  of  the  growth  delay  region,  the  model  behavior 
closely  describes  the  average  growth  rate  curves  (see  Figure  8). 

7.  BELOW  THE  SURFACE 

It  was  assumed  that  the  crack  growth  behavior  observed  on  the  surface 
was  Identical  to  that  whic..  was  occurring  within  the  specimen.  Failure 
of  the  specimen  revealed  that  following  the  initial  overload,  the  surface 
crack  movement  in  the  120  ksi  yield  strength  steel  specimen  lagged  the 
crack  movement  in  the  center  of  specimen  by  an  increasing  amount  with 
each  additional  overload  application.  Detailed  in  Table  IX  is  a  listing 
of  the  surface  and  Internal  crack  movements  as  related  to  applied  over¬ 
load  and  constant  amplitude  cycles  between  overloads.  The  prior  to  over¬ 
load  application  crack  position  on  the  surface  and  at  the  center  of  the 
s  ’ecimen  is  tabulated  In  Table  X. 
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TABU  VIII 


2 


Figure  8.  Model  Prediction  Applied  in  Reverse  Order  (Steady  State  then 
Nonsteady  state  rates)  Compared  to  Average  Behavior 
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TABLE  IX 

SURFACE  AND  INTERNAL  CRACK  MOVEMENTS  BETWEEN  SINGLE  OVERLOADS 


CRACK  MOVEMENTS  OBSERVED 

Overload 

£a** 

Mils 

H9I 

m 

KM 

Cycles  Between 
Overloads 

1A* 

35 

47 

2 

30,000 

1 

25 

39 

4 

35,000 

2 

29.5 

40 

3 

52,000 

3 

27.5 

38 

3 

42,000 

4 

23.5 

34 

2 

40,000 

5 

27 

34 

2 

42,000 

6 

40.5 

38 

2 

41,000 

7 

32 

38 

2 

46,000 

8 

26 

38 

3 

43,000 

AVG  (1  -  8) 

28.9 

37.4 

— — ...  — .  — .. 

2.5 

42,600 

*  Overload  1A  had  a  reduced  K*^  /K°°  of  1.8 

max  max 

**  an<^  were  the  observed  total  surface  movements, 

observed  total  internal  movements,  and  jumps  associated  with  the 
overloads. 
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Figure  9.  Schematic  of  No  Retardation  Boundary  Concept  to  Account  for 
Growth  Jumps 
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TABLE  X 

SURFACE  ANO  INTERNAL  CRACK  POSITION  PRIOR  TO  EACH  OVERLOAD  APPLICATION 


OVERLOAD 

CRACK  POSITION 

SURFACE 

In. 

INTERNAL 

In. 

1A* 

0.958 

0.958 

1 

0.993 

1.005 

2 

1.018 

1.044 

3 

1.047 

1.084 

4 

1.075 

1.122 

5 

1.099 

1.156 

6 

1.126 

1.190 

7 

1.166 

1.228 

8 

1.198 

1.266 

O.L.  No.  8  PLUS  43000  CYCLES 

1 .224 

1.304 

*  OVERLOAD  1A  HAD  A  REDUCED  K0L  /\f  of  1.8 

max  max 
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Note  is  made  in  Table  IX  of  the  incremental  jump  in  the  center  of 
specimen  crack  position  that  occurs  during  the  overload  cycle.  This 
jump  phenomenon  has  also  been  observed  by  Von  Euw  et  al  (Reference  4) 
and  by  Hall  (Reference  15).  The  average  crack  jump  during  the  overload 
was  approximately  2.5  mils,  approximately  two  orders  of  magnitude  higher 
than  what  would  Jiave  been  calculated  using  a  steady  state  growth  rate  at 

a  :.K  -  40  ksi  /in.  The  present  model  does  not  account  for  accelerated 
growth. 

It  is  suspected  that  a  generalization  of  the  use  of  the  load  inter¬ 
action  zone  concept  (Reference  3)  could  be  employed  as  a  first  step  in 
modeling  the  jump  phenomenon.  For  example,  consider  Figure  9  where  the 
base  line  no  retardation  curve  (K*)  has  been  established  for  a  low  level 
loading  which  produces  steady  state  growth.  A  high  level  load  is  applied 
at  crack  length  ag^,  establishing  a  subsequent  no  retardation  curve 
(«**).  Assume  that  the  incremental  movement  (accelerated  growth) 
accompanying  the  application  of  the  high  load  is  related  to  the  difference 
between  the  two  no  retardation  curves  (i.e.,  K*  and  K**).  Taking  the 
same  approach  that  was  employed  in  the  development  of  the  retardation 
model  presented  above,  one  could  set  the  residual  stress  intensity  factor 
equal  to  a  proportionality  factor  multiplied  by  the  difference  in  the 
no  retardation  stress  intensity  factors,  i.e. 

Kr  =  -A  (K**  -  K*)  (16) 

which  would  then  be  used  in  conjunction  with  Equation  1  to  calculate 
the  rate  of  crack  advance.  Insufficient  evidence  precludes  any  further 
discussion  of  the  applicability  of  the  load  interaction  concept  as 
applied  to  overload  induced  crack  growth  acceleration. 

Since  the  crack  growth  behavior  obtained  through  surface  observations 
is  significantly  different  from  that  noted  along  the  center  of  the 
failed  specimen,  It  Is  difficult  to  place  full  confidence  in  the  crack 
growth  rate  behavior  generated  using  surface  growth  data.  Possibly, 
there  is  a  rather  abrupt  deceleration  in  crack  advancement  following  an 
overload  but  this  question  will  have  to  be  addressed  by  future  investi¬ 
gations  of  the  delay  phenomenon. 
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Using  the  internal  crack  tip  position  (see  Table  X)  and  the  known 
applied  loads  employed  to  control  stress  intensity  factor  variables  on 
the  basis  of  surface  measurements,  it  was  possible  to  determine  the 
change  in  the  stress  intensity  factor  variables  associated  with  the  mid¬ 
specimen  position.  The  average  constant  amplitude  stress  intensity  range 
was  19  ksi  /in'.  It  changed  from  18  to  20  ksi  / in  between  the  first 

overload  and  the  last  constant  amplitude  cycle.  The  average  was 
_  r\i  r—— —  max 

21  ksi  J  in  ,  while  the  average  KUL  was  42  ksi  /  in.  Table  XI  provides 

max 

the  calculations  for  the  retardation  model  behavior  with  the  mid-stress 
intensity  factor  specimen  conditions.  Following  the-  numerical  integration 
through  the  overload  affected  zone  (a*  =  14.6  mils)  the  crack  was  allowed 
to  move  to  the  average  /.a position  (see  Table  IX).  The  calculated 
number  of  cycles  is  approximately  94  percent  of  the  average  number  of 
constant  amplitude  cycles  applied  between  overloads. 
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TABLE  XI 


RETARDATION  MODEL  CYCLE  CALCULATIONS  ASSUMING 
MID-SPECIMEN  STRESS  INTENSITY  FACTOR 


a' 

mils 

Aa 

mils 

— 

AN 

0 

1.95 

10.2 

8392 

1.95 

1.95 

11.5 

5717 

3.90 

1.95 

12.8 

4058 

5.85 

1.95 

14.1 

2978 

7.80 

1.95 

15.6 

2154 

9.75 

1.95 

17.2 

1577 

11.70 

1.95 

18.5 

1249 

13.65 

0.975 

19 

573 

14.62 

22.8 

19 

13406 

37.40* 

LAN  =  40104  *  Nm 

*  37.40  was  taken  as  average  movement  between  overloads  (see  Table  IX) 
which  could  be  associated  with  an  average  number  of  cycles  (N  )  of 

a  Vy 

42,600  to  produce  this  growth:  N^/NAVG  = 
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SECTION  IV 

SUMMARY 

Significant  differences  in  delay  cycles  were  observed  when  the  yield 
strength  of  a  4340  steel  was  varied  from  120  to  220  ksi.  Subjecting 
specimens  to  controlled  stress  intensity  factor  conditions:  single  peak 

overloads  with  =  40  ksi  /TrT followed  by  constant  amplitude  cycling 

w  max _ 

with  Kmax  =  20  ksi  /in  (R  =0.1),  indicated  that  for  the  same  external 
loading  {e.g.,  spectrum  loading)  more  retardation  can  be  expected  from 
a  steel  with  lower  yield  strength. 

Following  the  development  of  a  retardation  model  based  on  overload 
shut-off  and  threshold  stress  intensity  factor  conditions,  it  was  possible 
to  shew  that  a  lower  strength  steel  gives  more  retardation  for  a  given 
loading  condition  primarily  because  its  overload  created  load  interaction 
zone  is  significantly  larger. 

The  retardation  model  developed  herein  was  found  to  predict  to  with¬ 
in  10  percent  the  cycles  required  to  propagate  a  crack  from  its  position 
immediately  prior  to  overload  application  to  a  subsequent  position  a 
distance  one  plane  stress  plastic  zone  radius  ahead  of  its  pre-overload 
position.  Crack  growth  rates  obtained  from  surface  measurements  were 
not  adequately  predicted  by  the  retardation  model  but  note  was  made  of 
similarities  in  predicted  growth  rates  and  observed  behavior. 
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